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ABSTRACT 

A method for using the sub-millimeter band to determine photometric red- 
shifts of luminous high-z dusty galaxies is presented. It is based on the observa- 
tion that local ultra-luminous IR galaxies show an average spectral energy distri- 
bution which has a remarkably small dispersion for wavelengths A > 50/im. The 
shape of the long wavelength part is independent on whether the local ULIRG 
possesses an AGN or not. A local template is derived from a sample of ULIRG 
data presented by Klaas et al. (2001). Using this template and the internal 
dispersion in its shape, it is shown that observations at two fixed submm wave- 
lengths, e.g. 450/im and 850/xm, can be used to determine the photometric red- 
shift of galaxies in the redshift range 1 < z < 5. The uncertainty in the redshift 
arising from the template SED amounts to Az/z ~ 0.3, with an improvement at 
the higher redshift range. Implications for using this method together with the 
large instantaneous bandwidth of the Atacama Large Millimeter Array (ALMA) 
for deriving spectroscopic redshifts using CO emission is discussed. 

Subject headings: techniques: photometric, galaxies: distances and redshifts, 
galaxies: high redshift, galaxies: photometry, submillimeter 

1. Introduction 

The spectral shape of the far-infrared background suggests that as much as half of the 
energy ever emitted by stars and AGNs has been absorbed by dust grains and re-radiated 
in the far-infrared (Puget et al. 1996; Fixsen et al. 1998; Gispert, Lagache & Puget 2000). 
Recent observations at sub-millimeter wavelengths have revealed a population of objects 
believed to be dusty high redshift far-infrared luminous galaxies (Smail et al. 1997; Barger et 
al. 1998; Hughes et al. 1998; Eales et al. 1999). These sources emit most of their bolometric 
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luminosity in the far-infrared (FIR), and are believed to account for a large portion of the far- 
infrared background (e.g. Blain 1999). Although detailed redshift information is still lacking 
for most of these sources, they are believed to be situated at z > 1. Number-counts then 
indicate a large excess of sources compared to no-evolution models based on optical surveys 
(Guiderdoni et al. 1998; Blain 1999; Scott et al. 2002b), suggesting that the optically derived 
star formation density is significantly underestimated at high redshifts. This strong evolution 
is difficult to reproduce in semi-analytical hierarchical Cold Dark Matter models (Devriendt 
& Guiderdoni 2000; Kaviani, Haehnelt & Kauffmann 2002). Assuming similar dust and 
FIR properties between the high- and low-z objects, the CDM models under-estimate the 
observed number counts by almost two orders of magnitude (Kaviani et al. 2002). 

A remaining and central uncertainty in the discussion about the properties of the dusty 
high redshift population is their actual redshift distribution. Despite detection of >100 sub- 
millimeter sources, only a handful have secure optical or near-infrared counterparts (cf. Smail 
et al. 2002). This lack of counterparts shows that the submm sources are very faint (R>25 
and K>21), making it difficult to obtain redshifts even with 8-10m class telescopes, should 
the counterparts be identified. The presence of large amounts of dust also makes optical/NIR 
photometric redshift determinations unreliable, should the continuum be detected. 

Ultra-deep long wavelength radio continuum surveys have shown that many of the 
submm detected objects are related to weak radio sources (Smail et al. 2000). Chapman et 
al. (2002b) find that ~2/3 of sources detected at 850/im are also radio-identified at sub-mJy 
levels. For local galaxies there exists a very tight correlation between the long wavelength 
radio continuum and the far-infrared at 60/zm and 100//m, extending over almost four orders 
of magnitude in luminosity (Condon 1992). The physical mechanism for this correlation 
is believed to be the formation of massive stars, which heats the dust through their UV 
radiation, producing the FIR emission, and supernova explosions producing relativistic elec- 
trons which give synchrotron emission in the radio regime. Based on this tight correlation 
between non-thermal radio and thermal far-infrared fluxes for local galaxies, Carilli & Yun 
(1999) devised a technique to use the radio-FIR spectral index as a redshift indicator for 
distant submm sources. The technique uses the fact that for a fixed observed wavelength, 
the flux of synchrotron emission decreases with increasing redshift, while the opposite is true 
for the Rayleigh- Jeans part of the dust continuum. The first implementation used two local 
starburst galaxies as templates (see also Barger et al. 1998). A follow-up study (Carilli & 
Yun 2000) used 17 low-z star-forming galaxies as templates. The result is a submm-radio 
spectral index which depends on the redshift of the source. 

There are several caveats with this radio-FIR relation. The spectral index depends not 
only on the redshift of the source, but also on several other parameters, most notably the dust 
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temperature, and the radio continuum spectral index, for instance non-thermal vs. thermal 
contributions and the possible presence of an AGN. The decrease of the observed radio 
continuum with increasing redshift makes it increasingly demanding to detect the sources as 
z increases. The radio-FIR relation is therefore practical up to z ~ 2. Dunne et al. (2001) 
discuss the various dependencies, comparing with a local sample of FIR bright galaxies. 
Nevertheless, although the radio-FIR technique is too crude to be used as a photometric 
redshift indicator for individual sources, becomes unreliable at z > 2, and care has to be 
exercised in order not to include radio-loud AGNs, it has proven that the submm detected 
sources, seen as a class, are indeed situated at high redshifts. Very little is known about 
these sources from other wavelength bands, therefore this information has been crucial for 
establishing the presence of a high redshift population of dusty, luminous and, in the few 
known cases, massive galaxies. 

Several authors have speculated on the possibility of using two or more wavelength 
bands in the sub-millimeter as a mean of deriving photometric redshifts (e.g. Hughes et al. 
1998; Blain 1999; Blain, Barnard & Chapman 2002). The conclusion has been that this is 
unfeasable due to a degeneracy between the redshift and dust temperature (see Sect. 2.3). 
In a recent study by Hughes et al. (2002) the submm photometric method is revived using 
local starburst galaxies and a statistical treatment. 

In this paper an alternative method for determining the redshift of submm detected 
sources is presented, using only the Rayleigh- Jeans part of the dust SED. It is based on the 
existence of a template SED for local ULIRGs, which shows a remarkably uniform shape 
at wavelengths 60yum< A < 1250/im. Assuming that ultraluminous IR galaxies at high 
redshift have similar properties, the flux ratio of two submm bands, for instance at 850yum 
and 450/zm, is thus sensitive to the redshift of the source in the range 1 < z < 5. A 
comparison with submm flux ratios for a sample of high-z AGNs (mostly gravitationally 
lensed) shows that this technique is accurate enough to allow a redshift determination to 
within ~30%. With the advent of new submm telescopes and instruments, such as APEX 
(Atacama Pathfinder Experiment), ALMA (Atacama Large Millimeter Array) and SCUBA 
II, some of them capable of operating up to 325/im, the use of submm flux ratios is foreseen 
to become a technique to quickly estimate approximate redshifts of large numbers of submm 
detected sources. 



2. Dust Spectral Energy Distributions 

In this section I will briefly discuss dust spectral energy distributions and how these 
can be used to describe physical properties of unresolved dust sources observed in the IR- 
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to-submm. The main conclusion is that extreme care has to be exercised when doing this 
interpretation since it can easily lead to wrong inferences about dust mass and the physical 
size of the emitting region. 



2.1. Models 

The far-infrared spectral energy distribution of dust emission is purely thermal and 
usually represented by a modified blackbody curve, S v oc v^B^Ta), where B v is the Planck 
function, T d the dust temperature and u 13 is the grain emissivity power-law index. The index 
is in the range /3 — 1 — 2. A single-size and uniform grain composition would lead to f3 — 1. 
When /3 > 1 it is indicative of a distribution of grain sizes and/or compositions. Such simple 
representation has successfully been used for cold dust components where a large part of 
the SED is optically thin. When r ps 1 or larger, the observed dust emission needs to be 
described by an expression of the form: 

S v = Q s B u (T d ) (l - e- T ») , (1) 

where Q s is the solid angle of the source emissivity distribution and r„ is the opacity of the 
dust. Setting t v = {yjv^f gives S u oc v^B v i (T d ) for r v <C 1 and S v oc B v {T d ) for r v ^> 1. 
The critical frequency uq is the frequency where t v — 1. The FIR luminosity is obtained by 
integrating Eq. 1 over all frequencies: 

oo 

L m = A7r(l + z) 3 D 2 A Js u du , (2) 

o 

where Da is the angular size distance 1 With the substitution x = hu/kT, we can express 
Eq. 2 in a form directly accessible for integration (e.g. Wiklind & Alloin 2003): 



L© 



fx 3 (l- e~( ax ^) 

= 8.53 x 10 10 (1 + zf D\T* tt s / — i J -dx, 

J e x — 1 



(3) 



The integral can be evaluated numerically with appropriate values of the parameter a = 
kT d /hu . The solid angle Q s is a parameter derived in the fitting procedure. In the event 
of a single dust component and a source unresolved by the telescope, Q s can be estimated 
from the measured flux S Vl at a given restframe frequency v r 



B Vr (T d ) (l-rW^) 



1 Unless otherwise stated, all distances will be expressed in Mpc, frequencies in GHz and fluxes in Jansky. 
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« 6.782 x lO -4 ^ -3 ^ 

A rough estimate of the equivalent spherical radius of the emitting region can be obtained 
from Q s ; r D^^jQ s /n. 

An estimate of the dust mass from the infrared flux requires either optically thin emis- 
sion combined with a knowledge of the grain properties, or optically thick emission and a 
knowledge of the geometry of the emission region (cf. Hildebrand 1983). The grain proper- 
ties are characterized through the macroscopic mass absorption coefficient, k v (cf. Hughes, 
Dunlop & Rawlins 1997). Combining several estimates (cf, Hildebrand 1983; Hughes et al. 
1997), the mass absorption coefficient can be described as: 

^°- 15 (ra) 1 ' 5 ^kg" 1 • (5) 

This expression corresponds to a grain composition similar to that found in the Milky Way. 
At frequencies where the emission is optically thin, the dust mass can now be determined 
from 

M d = - l u ° b : T M i + Z )\ (6) 

which can be expressed as: 

/ y \ -9/2 

M d ~ 4.08 x 10 4 S Uobs D\ (^^J {e hUr/kTi - l) (1 + zf M . (7) 



1 — e -(^r/^o) /3 



(4) 



2.2. Multiple components 

In reality, it is likely that a galaxy contains several different dust components each 
characterized by a different set of parameters. Lacking observational data with sufficient 
spatial resolution, it is usually assumed that the dust emission can be decomposed into 
three main components (e.g. Rowan- Robinson 1986, 1992); cool dust associated with a 
cirrus component heated by the diffuse interstellar radiation field, warm dust associated 
with regions of massive star formation, and, if the galaxy harbors an AGN, a hot component 
representing dust heated by the central AGN. These three dust components will dominate 
the submm, FIR and mid-IR wavelength regimes, respectively. There are, however, no 
assumptions about their actual sizes or their geometries. 

More elaborate models solve the radiation of transfer in a dusty medium containing 
heating sources (e.g. Granato, Danese & Franceschini 1996; Efstathiou, Rowan-Robinson & 
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Siebenmorgen 2000). Due to the large number of parameters, involving dust variables as a 
function of the geometry, these methods can only be safely applied to sources and parts of 
the SED which are known to be associated with a single heating source, e.g. an AGN. 

The representation of the dust SED discussed in Sect. 2.1 involves four parameters per 
component. Since most observed SEDs contain only a few data points, a fit of Eq. 1 is prone 
to a relatively large degree of freedom, especially if more than one dust component is used. 
Moreover, at long wavelengths there exists a degeneracy between the dust temperature T d 
and the grain emissivity power-law index (3 (e.g. Blain et al. 2002). This degeneracy is 
difficult to resolve without an observationally well sampled SED. 

Due to the observationally sparsely sampled dust SEDs, the dust emission is often char- 
acterized by fitting a single component. This usually gives a reasonable fit to the observed 
data points, but should not be over-interpreted in terms of physical conditions of the dust 
grains or of global dust parameters, especially dust mass. This is illustrated in Fig. 1 where 
the dust SED of the local ULIRG 17208-0014 has been fitted using both one and two dust 
component SEDs of the type described above (see also Klaas et al. 2001). The single com- 
ponent SED is characterized by a dust temperature of 59 K and an exponent f3 = 1.8. The 
two component SED has warm dust grains characterized by a temperature of 53 K and 
P = 1.6, and a cool component with T d = 27 K and (3 = 2.0. While the two component 
SED has critical frequencies z/ corresponding to 35/zm and 51/im, respectively, the single 
component SED becomes optically thick already at 214/im. Both models give approximately 
the same L F ir, showing that they describe the shape of the SED reasonably well. The sin- 
gle component SED, however, gives a total dust mass which is 3.5 times smaller than the 
two-component SED. This is due to the highly non-linear relation between Md us t/LFiR as 
a function of dust temperature; while Md oc exp (hu r /kTd), (approximately oc T^ 1 at long 
wavelengths), the luminosity L F ir oc T^ -5 . It is thus possible to 'hide' a large amount of 
dust in a low temperature component. In terms of typical size scales, the cool dust extends 
over an area which is ~150 times larger than the warm component. 

Although the resulting x 2 values for the fits shown in Fig. 1 are similar, they describe 
dust components with widely different physical properties. The two fits, however, both give 
a satisfactory description of the overall shape of the SED. This is important to keep in mind 
when discussing and comparing dust properties of different types of objects, and, as argued 
here, when constructing a template SED. 
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2.3. Degeneracy between T d and redshift 

Lowering the temperature of a modified blackbody curve is equivalent to shifting the 
entire spectrum towards longer wavelengths. A similar shift can be achieved by keeping the 
temperature constant but moving the source to a larger redshift. Hence, without knowing the 
dust temperature apriori, it is impossible to distinguish between a source at large redshift 
with a given Td and a source at low redshift with a dust SED characterized by a lower 
temperature (e.g. Blain 1999; Eales et al. 1999; Blain, Barnard & Chapman 2002). 

In order to use two flux densities at the FIR/submm part of a dust SED to determine 
a photometric redshift of the source, it is thus necessary to know the shape of the SED in 
sufficient detail. This can be done if the global physical parameters for the emission can 
be established. If this is not the case, the degeneracy between the dust temperature and 
redshift makes progress essentially impossible. 

The effect of this degeneracy is illustrated in Fig. 2 for a simple dust SED. The ratio of 
flux densities at 325/im and 850/im of a single component modified blackbody function (given 
by Eq 1) are plotted as a function of redshift. The result for six different dust temperatures 
between 30 K and 80 K are shown. Despite the fact that the gradient of the flux ratio as a 
function of redshift is quite steep, the effect of the dust temperature precludes its use as a 
redshift indicator unless T d is known. 

Attempts to characterize local FIR luminous galaxies in terms of dust parameters show 
a relatively broad range in T d and (3 (e.g. Sanders & Mirabel 1996; Klaas et al. 2001). 
As shown above, a sparsely sampled dust SED can be equally well modeled with one and 
two component dust SEDs, with different sets of dust parameters, leading to an even larger 
spread in parameter values for different sources. Similarly good fits can also be achieved for 
different dust temperatures by varying the index f3. The shape, however, remains basically 
identical and the ratio of FIR/submm flux densities are approximately the same for the 
various fits. 



3. Photometric Redshift 

3.1. The local template 

Local ULIRGs are characterized by a dominating FIR luminosity, LpiR > 10 12 L , large 
dust masses, large molecular gas masses and signs of strong gravitational interaction and/or 
merging (e.g. Sanders & Mirabel 1996; Rigopoulou et al. 1999). In addition, the dust and 
gas is strongly concentrated to the central regions of the galaxies (Downes & Solomon 1998). 
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The heating source of the large FIR luminosity is in many cases pure star formation, but 
a significant fraction of all ULIRGs also contain an AGN (Rigopoulou et al. 1999; Klaas 
et al. 2001). Local ULIRGs have properties which in many respects are similar to those 
of the submm detected population, especially those few that have known redshift and have 
been studied in greater detail (cf. Ivison 2000). Two submm detected galaxies are known 
to contain large amount of molecular gas as seen through their CO emission (Frayer et al. 
1998, 1999). Using the standard conversion ratio between CO and molecular hydrogen, H 2 , 
the gas masses are in the range (0.5 — 2.0) x 10 11 M , showing that these objects are indeed 
as massive and gas-rich as local ULIRGs. Local ULIRGs are therefore the most likely low-z 
counterparts to the high redshift dusty submm objects. 

A recent compilation of NIR-to-submm observations of 41 local ultraluminous infrared 
galaxies was presented in Klaas et al. (2001). The sample consists of galaxies with L FIR > 
10 12 L , with most objects having redshift z < 0.1. The observational data was obtained 
with several different instruments. The NIR bands (1.2 and 2.2/zm) was observed using 
the 2.2m Calar Alto telescope. The IR data between 10-200/im was obtained using the 
ISO. Sub-millimeter data at 450/im and 850/zm was obtained for a sub-set of 16 objects 
using the JCMT. A few millimeter continuum observations were also done with the SEST 
(Swedish-ESO Submillimeter Telescope). The SEDs represent the most complete set of 
IR/FIR/submm photometric data obtained for local ULIRGs to date. Given the observed 
similarities between submm detected sources and the local ULIRGs, the Klaas et al. sample 
is well suited for creating a 'template' SED to be used for the high redshift counterparts. 

In October 2002 we observed 8 of the sources in the Klaas et al. (2001) sample using 
the SIMBA bolometer array at the Swedish-ESO-Submillimeter telescope (SEST) in Chile. 
The SIMBA instrument is a 37-channel bolometer array operating at 1250/im (see Nyman 
et al. 2001 for a description of SIMBA). It uses a fast mapping technique and achieves a 
sensitivity of ~70mJys -1 / 2 pixel -1 . Four of the local ULIRGs were positively detected. The 
results are presented in Table 1. The 1250/im result for the source 17208-0014 is used in 
Fig. 1, where the long wavelength data implies the presence of a cold dust component. This 
aspect of the dust continuum will be discussed in a forthcoming paper. The 1250/im data is 
not used in the present analysis. 

In order to avoid redshift corrections to the local sample, 4 objects from the Klaas et al. 
(2001) source list, with redshifts in the range z = 0.2 — 0.3, were discarded. The remaining 
37 objects all have z < 0.1. The SEDs of these 37 sources, as tabulated by Klaas et al., are 
plotted in Fig. 3. No corrections or normalization have been done to the SEDs at this stage. 
It is clear from the figure that although the observed flux densities vary by two orders of 
magnitude, the general shapes of the SEDs appear quite uniform. 
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In order to quantify the impression of a uniform shape of the SED, the individual SEDs, 
Si where multiplied with a constant Ai chosen such that the overall dispersion is minimized. 
In the present case, all the wavelengths were given equal weights. No significant change could 
be found when restricting the wavelengths to the FIR and sub-millimeter. The result is shown 
in Fig. 4, where the average flux (arbitrary normalization) at each observed wavelength is 
shown together with the la deviation. In Fig. 4 the resulting individual measurements are 
also shown. It is clear that the impression of a general shape of the SEDs is correct only for 
wavelengths A > 50yum. At mid-IR wavelengths, the dispersion is considerable. At near-IR 
wavelengths the dispersion is again smaller than in the mid-IR. 



3.2. The method 

3.2.1. Fitting a single SED 

The small dispersion of the average spectral energy distribution of local ULIRGs at wave- 
lengths longer than 50/zm (Fig. 4) allows a modified blackbody curve of the type discussed 
in Sect. 2 to be fitted to the average SED over this interval. Although a two component 
fit is marginally superior to a single component, especially around 90 — 120/im, a single 
component dust SED was used for simplicity. The fitted function is given by Eq. 1 and the 
uncertainty associated with each observed wavelength is taken to be the dispersion obtained 
from the the average SED shown in Fig. 4. 

A least-square fit of Eq. 1 gives the following parameters: Q = (1.6 ± 0.6) x 10 -13 , 
P = 1.8 ± 0.5, u = (1.2 ± 0.4) x 10 12 Hz (250/W) and T d = 68 ± 9 K. The Q parameter 
is necessary for the fit but otherwise arbitrary since the average SED has already been 
normalized. The x 2 is 0-89. As emphasized in Sect. 2, the physical reality of the dust 
properties associated with the fit is debatable. However, the shape of the SED is well 
described by the fitted parameter. 

The resulting fitted SED is shown in Fig. 5 together with the normalized data points. 
Only data in the wavelength range 60 < A < 850/im were used in the fit. The fitted curve 
coincides very tightly with all the average data points in the range 25 < A < 850/im, but 
only the 60 — 850/im range will be considered in the following. The 1250/im data were 
omitted because of a small number of observations and the possibility of contribution from 
a very cold dust component. Such a cold component would not be visible when observing 
at shorter wavelengths and definitely not for sources at high redshifts. The presence of such 
a cold dust component in at least one of the sources in the present sample is suggested by 
the excess flux density at A = 1250/xm seen in Fig. 1. 
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3.2.2. Monte Carlo simulation 

In order to estimate the uncertainty associated with the fit of the single dust component, 
given the dispersion of the average SED, a Monte Carlo simulation was done. For 10 4 real- 
izations the average values of the underlying data points were stochastically and uniformly 
varied within ±lo\ For each simulated data set a least-square fit of Eq. 1 was done. The fits 
were done using data in the wavelength range 60/im < A < 850/im. The resulting SEDs are 
shown in Fig. 6. All the 10 4 curves are plotted and they define an envelope which remains 
well constrained in the wavelength range 60 — 850/im, but which diverges at shorter and 
longer wavelengths (where the fits are no longer constrained by data). The fit to each of the 
Monte Carlo realizations where done in the same manner as the fit discussed above. 



3.2.3. Flux ratios 

The observed flux at a given wavelength band is obtained by multiplying Eq. 1, using 
the fitted parameters, with a normalized filter response function and then integrating over 
wavelength and normalizing with the filter response function. Here a simple top-hat function 
is used, with the same width as the filters used on SCUBA, i.e. 30 GHz for both 450/im and 
850/im. The same width was also used for the 325/im band 2 . By redshifting the underlying 
SED and correcting for the (1 + z) compression, the flux ratio at a given redshift was 
calculated. 

The flux density ratio at two fixed wavelengths as a function of redshift was obtained 
for each of the Monte Carlo realizations. They are plotted in Fig. 7 for both 325/850/im and 
450/850/im. The grey-scale corresponds to the density of curves passing through a given 
point in the diagram, the darker the shade the higher the density. The white full drawn 
line marks the ridge of highest density and it is flanked by two lines marking the boundary 
containing 95% of the curves on either side of the ridge. By taking a cut at a constant flux 
ratio, it is possible to assess the formal error in estimated redshift. This error is for the most 
part non-gaussian. In reality there will be an uncertainty associated with the measured flux 
ratio as well. The predicted redshift will be enclosed by a box-shaped region, limited by 
the high end of the flux ratio uncertainty at the lower redshift end and the lower flux ratio 
uncertainty at the higher redshift end. 



2 This wavelength was chosen in anticipation of the planned bolometer camera for APEX (Atacama 
Pathfinder Experiment). The precise wavelength coverage and the shape of the filter response function arc 
not known presently. The results presented here for this short wavelength band must therefore be recalculated 
once the filter response is determined. 
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At redshifts z < 1, the dispersion in redshift is too large for the flux ratio to contain any 
predictive power of the redshift. This is due to the fact that at small redshifts, both of the 
observed wavelengths probe a section of the Rayleigh- Jeans part of the modified blackbody 
curve with a close to a pure power law form. The flux ratio therefore remains almost constant. 
The insensitivity to low redshifts is less pronounced for the 325/850/im flux ratio since the 
shorter wavelength will be affected by non-linearities in the dust SED at lower redshifts than 
the 450/xm band. At high enough redshifts the shorter wavelength band will reach the Wien 
side of the modified blackbody curve. The flux ratios will again start to become more or less 
constant. For the present sample this effect will set in at redshifts z > 5. The most accurate 
predictive power for the redshift therefore exists in the redshift range 1 < z < 5. 



4. Discussion and Conclusion 

4.1. Selection effects 

Any method of photometric redshift is based on the assumption that we know or can 
model the spectral energy distribution of high redshift objects. Modeling the FIR/submm 
SED is difficult due to a degeneracy between the dust temperature and the grain emissivity 
power-law index (3. Moreover, as discussed in Sect. 2.3, lowering the dust temperature of a 
nearby source gives a similar effect on the SED as moving the source to a higher redshift 
while keeping the temperature constant. These degeneracies makes it impossible to use the 
FIR/submm SED for photometric redshift determinations if one do not have a sufficiently 
good apriori knowledge of the dust properties of the object in question. There appears to 
exist a relation between the dust temperature and FIR luminosity (e.g. Dunne et al. 2000; 
Blain et al. 2002), which can possibly be used to break the degeneracy. This Luminosity- 
Temperature (LT) relation is, however, presently poorly determined. 

In Sect. 3.1 it was shown that the SED of local ultraluminous IR galaxies has a remark- 
ably similar shape at wavelengths longer than A > 50/im. The similarity in the shape of 
the SEDs could very well be an effect of sample selection. The local galaxies was selected 
as being ULIRGs, meaning that they have a Lj?ir > 10 12 L . At the same time, all the 
sources have Lpm < 10 13 L . If there exists a LT- relation, this selection of sources within a 
limited luminosity range favors those which have similar dust temperatures. Whenever de- 
tailed properties are known for the high redshift submm detected sources, which is the case 
for 2-3 objects, they have properties very similar to the local ULIRGs in terms of Lfir, gas 
mass and star formation. If the submm detected objects are at high redshift they must have 
L F ir exceeding 10 12 L , or they would not be detected with present day instrumentation. 
They are therefore likely to occupy the same region of the LT-relation as the local ULIRGs. 
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Making the assumption that the high redshift population as a whole is similar to the local 
ULIRGs, the local sample can be used as a template. 

While the local sample shows a remarkably small dispersion in the shape of the long 
wavelength part of their SEDs, they differ considerably at mid-IR wavelengths. In this 
wavelength regime the IR emission is likely to show optical depth effects. The resulting 
SED is then dependent on the geometry of the emitting region and its orientation relative 
to the observer. The presence of an AGN would also be most apparent in this wavelength 
regime (see Klaas et al. 2001). While the shape of the SED at mid-IR contains substantial 
diagnostic power for the presence of a compact heating source for the dust, whether it is 
massive star formation or an AGN, it does not affect the use of submm bands for photometric 
redshifts as long as z < 10. 

4.2. Comparison with other methods 

The assumption made here of the existence of a template SED is exactly the same as 
done when using the radio-FIR spectral index method (Carilli & Yun 1999, 2000). The main 
difference between the two method is that the FIR SED method does not need to exclude 
radio- loud AGNs, and it works out to redshifts z ~ 5. This is important since it is very 
difficult or impossible to distinguish between high redshift low-luminosity AGNs and non- 
AGNs among /iJy radio continuum sources without milli-arcsecond imaging. One example 
of a high redshift radio-loud AGN, detected both as a submm source and in radio continuum 
is the gravitationally lensed H1413+117. The observed submm flux is much lower than 
expected from the radio continuum flux (Carilli & Yun 1999). The sample used here consists 
of a mixture of pure starburst and AGNs systems (Klaas et al. 2001); 8 are designated as 
starburst systems, 13 as LINERs, 8 as Seyferts (2 Sy 1 and 6 Sy 2), and 8 remain unassigned 
of a spectral class. As discussed above, the presence of an AGN is most evident in the mid-IR 
range and not at A > 60/im. 

An alternative photometric redshift method for use with FIR and submm observations 
was presented by Hughes et al. (2002). This method use a large number of template SEDs, 
containing both starburst systems and AGNs, spanning a large range of intrinsic luminosities. 
The local IRAS 60/im luminosity function is evolved with redshift to explain the observed 
number-counts. The luminosity function is then randomly populated by the SEDs from the 
template catalog. Photometric redshifts are determined by calculating the probability that 
the colors of an observed submm source are consistent with the colors of every template SED 
at a given redshift. No correlation between the SED, luminosity and redshift is assumed. 
While models of this kind incorporate lower luminosity objects, the existing instrumentation 
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limits detection of high redshift objects to the most luminous ones. Until submm observations 
can routinely be done with rms noise levels of 0.1 mJy or lower, and in view of a possible 
luminosity-temperature relation, it seems prudent to restrict the template to objects with 
luminosities large enough to be detectable at high redshifts. 

4.3. Comparison with observations 

The flux density ratio 450 / um/850 / u.m, derived from the template, as a function of redshift 
is shown in Fig. 8. Also shown are measured flux density ratios for high-z sources with known 
spectroscopic redshifts (see Table 2). Six of these are taken from a survey of gravitationally 
lensed AGNs (Barvainis & Ivison 2002). Three are submm detected galaxies with known 
redshift: SMM02399-0136 at z = 2.81 (Ivison et al. 1998), SMM02399-0134 at z = 1.06 and 
SMM1411+0252 at z = 2.56 (Ivison et al. 2000). Two sources detected by ISO at 170/im 
(Chapman et al. 2002): FN1-40 at z = 0.45 and FN1-64 at z = 0.91. Two additional 
gravitational lenses are also included, F10214+4724 at z — 2.3 (Rowan-Robinson et al. 
1993) and APM08279+5255 at z = 3.91 (Irwin et al. 1998). The ISO sources are situated 
at relatively small redshifts and still have large uncertainties associated with their observed 
450/im flux densities. Their predicted redshifts with the present photometric method are 
therefore uncertain. Nevertheless, the present method achieves a better result than the 
radio-FIR method which predicts redshifts of 1.1 and 1.3 instead of the spectroscopically 
derived 0.45 and 0.9. 

Most of the sources shown in Fig. 8 and listed in Table 2 are gravitationally lensed 
and it is possible that differential magnification can distort the flux density ratios. This can 
happen if the emission regions dominating at 450/im has a smaller spatial extent than that of 
850/iin. The effect would be to enhance the 450/im flux relative to that at 850/im and move 
the data points upwards in the diagram. This, however, is only expected to be of concern 
for the most strongly lensed sources. Indeed, the data points deviating significantly from 
the expected ratios are due to the strongest lenses (F10214+4724 and APM08279+5255). 
Another deviant source is SMMJ02399-0134 with a spectroscopic redshift of z — 1.06, where 
the indicated photometric redshift is 2.2t{j;g (see Table 2). This, however, is a redshift where 
the present method is less accurate due to the relatively small redshift. Apart from these 
three sources, the agreement with the redshift expected from the template model is good. 

The low redshift sample of Dunne et al. (2001) generally fall below the predicted flux 
density ratio (given their redshift). This sample, however, only contains two sources which 
can be classified as ULIRG and hence, is not expected to conform to the ULIRG template. 
In particular, their average dust temperature is considerably lower (Dunne et al. 2001), in 
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agreement with a possible luminosity-temperature relation. 

The spectroscopic and photometric redshifts are given in Table 2. The photometric 
redshifts are given with the errors calculated as a convolution of the photometric uncer- 
tainty and the 95% confidence limit derived from the template SED. The spectroscopic and 
photometric redshifts are compared graphically in Fig. 9. The uncertainty is dominated by 
the errors associated with the 450/zm flux densities rather than the spread in template flux 
ratios. 



4.4. Photometric to spectroscopic redshifts 

Ultimately it is preferable to determine the redshift of all the submm detected galaxies 
spectroscopically. As mentioned in Sect. 1 this is in most cases very difficult due to the 
very weak emission at optical and NIR wavelengths. An alternative method will become 
available when ALMA (Atacama Large Millimeter Array) comes on-line. The sensitivity of 
this interferometer will be high enough to allow detection of CO emission from high redshift 
objects in a relatively short observing time (Blain et al. 2000). The instantaneous wavelength 
coverage will be much larger than present day instruments. The projected bandwidth is 
16 GHz, corresponding to a redshift intervall Az/z 0.2 — 0.3 in the 3mm atmospheric 
window. Although still relatively limited, it represents an order of magnitude improvement 
to existing millimeter and submillimeter instruments. The 3mm atmospheric window is ideal 
for searching for high redshift CO emission. Except for a redshift intervall 0.44 < z < 1.0, 
there will always be at least one CO rotational transition within 80 — 115 GHz. For redshift 
z > 3 there will be at least two transitions within the band. 

If the photometric errors can be reduced to such a level that the uncertainty in the 
photometric redshift is dominated by the dispersion in template SED, the redshift uncertainty 
would be ~30% for redshifts z > 1.5. This is illustrated in Fig. 10, where Az/z is plotted 
as a function of redshift for the 95% and 90% confidence limit of the template SED. Also 
shown in Fig. 10 are the corresponding Az/z values of four CO transitions as observed in 
the 3mm band with the projected ALMA receiver bandwidth of 16 GHz. The 3mm band 
used here stretches from 80 GHz to 116 GHz, although the first receiver implementation 
is likely to cover a somewhat smaller band. Using higher frequencies results in a smaller 
redshift coverage and does not cover the rotational transitions believed to give the strongest 
signal (e.g. Combes, Maoli & Omont 1999). For the highest redshifts covered by the 3mm 
band, and thus the higher J-transitions, the redshift uncertainty in the photometric method 
is comparable or smaller than the redshift coverage of the ALMA bandwidth. For lower 
redshifts it will be necessary to use two tunings in order to brackett a photometrically 
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estimated redshift. Improvements of the photometric redshift method may remedy this, 
making it possible to use a single tuning of ALMA when searching for redshifted CO emission 
even for lower redshift sources. 



5. Conclusion 

Despite a serious degeneracy between the dust temperature and the redshift for a mod- 
ified blackbody curve, it is shown that an average dust spectral energy distribution, based 
on a local sample of ultraluminous IR galaxies, can be used for photometric redshift deter- 
mination in the intervall 1 < z < 5. 

The local sample of ultraluminous IR galaxies (Klaas et al. 2001) shows a remarkably 
small dispersion in the shape of the SED at wavelengths longward of A = 50/im. The local 
sample contains both starburst and AGN dominated sources. 

Using this template and the internal dispersion in its shape, it is shown that observations 
at two fixed submm wavelengths, e.g. 450/zm and 850/zm, can be used to determine the 
photometric redshift of galaxies in the redshift range 1 < z < 5. The uncertainty in the 
redshift arising from the template SED amounts to Az/z ~ 0.3, with an improvement at the 
highest redshift range. 
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Table 1. SIMBA data at 1250/iin 



Name 


Sl250 

mJy 


a 
mJy 


Comments 


00199-7426 


49 


13 




05189-2524 


11 


3 




06035-7102 


<43 




3a 


06206-6315 


<13 




3a 


17208-0014 


64 


13 




20100-4156 


21 


6 




E148IG002 


<16 




3a 


23230-6726 


<21 




3a 



Note. - Sources from the Klaas et al. 
(2001) sample observed with the SIMBA 
bolometer array at the SEST in October 
2002. 
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Table 2. Submm flux data 



Name 


z spcc 


/850 


°"850 


/450 


C450 




R a 




z phot 


Reference 








(mJy) 


(mJy) 


(mJy) 


(mJy) 












SMM J14011+0252 


o 
Z 


.00 


12.3 


1.7 


42 


7 


3. 


.4 ± 0. 


.7 


2 4+2-1 

— 0.5 


(1) (2) 


SMM J02399-0136 


o 
Z 


.01 


23.0 


1.9 


85 


10 


3. 


.7±0. 


.5 


9 9+0.4 
z - z -0.3 


(1) (2) 


HR10 


1 
1 


A A 


8.7 


1.6 


<180 






<21 




- 


(3) 


SMM J02399-0134 


1 
1 


.UD 


11.0 


1.9 


42 


10 


3. 


.8±1. 


.1 


9 9+O.9 
z - z -0.6 


(1) 


FN1-64 


n 
u 


Q1 

.y i 


5.9 


1.4 


50 


19 


8. 


.5 ±4. 


.0 


0.2t L6 


(4) 


FN1-40 


A 

u 


/I tz 
.40 


6.3 


1.4 


42 


30 


6. 


.7 ±5. 


.0 


0.8+ 37 


(4) 


SMM J09429+4658 






17.2 


1.9 


61 


7 


3 


.5±0. 


.6 


9 0+0.4 
z.o_ 4 


(1) 


SMM J14009+0252 






15.6 


1.9 


33 


9 


2. 


.1 ±0. 


.6 




(1) 


SMM J00266+1708 






18.6 


2.4 


<60 






<3.2 




> 2.6 


(5) 


Lockman850.1 






10.5 


1.6 


35 


10 


3. 


.3 ± 1. 


.1 


9 c+1.1 


(6) 


tiub 850.1 






7.0 


0.4 


<zl 






<3.0 




> z.o 


{') 


2237+0305 


1 


.70 


2.8 


0.9 


<17 






<6.1 




> 1.0 


(8) 


HE0230-2130 


2 


.16 


21.0 


1.7 


77 


13 


3 


.7±0. 


.7 


9 9+O.5 
z.Z_ 4 


(8) 


F10214+4724 


2 


.29 


50.0 


5.0 


273 


45 


5. 


.5±1. 


.1 


I Q+0.5 

■■■• o -0.4 


(8) 


H1413+117 


2 


.56 


58.8 


8.1 


224 


38 


3. 


.8±0. 


.8 


9 9+O.6 
z - z -0.5 


(8) 


MG J0414+0534 


2 


.64 


25.3 


1.8 


66 


16 


2. 


.6±0. 


.7 


q 9+I.O 
J - z -0.6 


(8) 


UM673 


2 


.73 


12.0 


2.2 


<40 






<3.3 




> 2.5 


(8) 


RX J0911.4+0551 


2 


.81 


26.7 


1.4 


65 


19 


2. 


.4±0. 


.7 


3.4+o:? 


(8) 


MG J0751+2716 


3 


.21 


25.8 


1.3 


71 


15 


2. 


.8±0. 


.6 


3.ol8:I 


(8) 


CLASS B1359+154 


3 


.24 


11.5 


1.9 


39 


10 


3. 


.4±1. 


.0 


2 4 +L0 

Z ' 4 -0.7 


(8) 


APM08279+5255 


3 


.91 


84.0 


3.0 


285 


11 


3. 


.4±0. 


.2 


2 4 +0 - 2 


(8) 



a R = S'(450^m)/S'(850^m) 

Note. — References: (1) Smail et al. (2002) (2) Ivison et al. (2000); (3) Cimatti et al. (1998); (4) 
Chapman et al. (2002); (5) Frayer et al. (2000); (6) Lutz et al. (2001); (7) Hughes et al. (1998); (8) 
Barvainis & Ivison (2002) 
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Fig. 1. — a) Fit of a one-component modified balckbody function to the observed flux 
densities of the local ULIRG 17208-0014. The parameters for the fit are: Q = (9.1 ± 1.1) x 
1(T 14 , (3 = 1.8 ± 0.2, u = (1.4 ± 0.2) x 10 12 Hz, T d = 59 ± 1.4 K. b) Fit of a two-component 
modified blackbody function to the same flux density values. The fit is as good as the one- 
component and results in the following values: Q = 34.4 x 10~ 14 , (3 — 1.6, i/ — 8.6 x 10 12 
Hz, T d = 53 K, for the warm component, and Q = 371 x 10~ 14 , (3 = 2.0, u = 5.9 x 10 12 Hz, 
T d = 27 K, for the cold component. 
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Fig. 2. — The flux ratio at fixed observed wavlengths 325/!m/850/im for a modified blackbody 
spectral energy distribution as a function of redshift of the source. The flux ratios for 6 
different dust temperatures are shown, illustrating the degeneracy between redshift and dust 
temperature for a given flux density ratio. The dust SED is represented by: f u oc B u (T d ) : 
where (3 = 1.8. 
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Fig. 3. — The spectral energy distributions of the 37 ULIRGs part of the local sample (from 
Klaas et al. 2001). The observed fluxes are plotted without any corrections or normalization. 
Despite a range of fluxes over two orders of magnitude, the general shape of the SEDs remains 
the same. 
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Fig. 4. — The spectral energy distributions of the 37 ULIRGs part of the local sample (from 
Klaas et al. 2001). Here the individual SEDs have been multiplied by a constant factor, 
chosen such that the overall dispersion is minimized. The individual data points, the average 
as well as the la dispersions are shown (light grey). The line just connects the average flux 
measurement at each observed wavelength. 
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Fig. 5. — A fit of a modified blackbody curve (see Eq. 1) made to the average SED obtained 
by minimizing the overall dispersion. The fit has been made over the wavelength intervall 
60-850/im. The best-fit parameters are: Q = (1.6 ± 0.6) x 10 13 , (3 = 1.8 ± 0.5, u Q = 
(1.2 ± 0.4) x 10 12 Hz and T d = 68 ± 9 K. 
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Fig. 6. — Fits of modified blackbody curves resulting from letting the flux density values at 
each wavelength vary uniformly ±1(7. The grey area represents 10 4 simulated fits. 
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Fig. 7. — The ratio of two fixed wavelength bands, each of width 30 GHz, as a function of 
redshift of the emitting source. The grey area corresponds to the results from 10 4 Monte 
Carlo simulations and represents the density of curves passing through a given point (dark 
means high density). The maximum density is outlined by the white curve and it is flanked 
by two curves representing the boundary containing ±95% of the curves, a) The ratio of 
325 / um/850/im and b) 450/mi/850/mi. The 450yum and 850/iin bands are chosen to represent 
the existing SCUBA filters. 
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Fig. 8. — The flux ratio of 450/im/850/im, measured at fixed wavelenghts, versus redshift of 
the emitting source. The results from 10 4 Monte Carlo simulations are shown. The greyscale 
represents the density of curves passing through a given point (dark means high density). The 
maximum density is outlined by the white curve and it is flanked by two curves representing 
the boundary containing ±95% of the curves. Also shown are the existing sources observed 
at both 450/iin and 850/xm and with known redshift. The open triangles are taken from 
a sample of high redshift QSOs (Barvainis & Ivison 2002). All of these are gravitationally 
lensed but differential magnification is likely to be significant only for the strongest lenses 
(IRAS F10214+4724 and APM08279+5255). Filled circles correspond to submm detected 
galaxies. 
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Spectroscopic z 

Fig. 9. — Comparison between photometric redshifts using the present method and spec- 
troscopic redshift for a sample of submm detected objects and high-z quasars (see Table 2). 
The full drawn line corresponds to a one-to-one relation between photometric and spectro- 
scopic redshifts and is not a fit to the data. The errors in photometric redshifts have been 
derived using the observed photometric errors convolved with the 95% confidence level for 
the template. 
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Fig. 10. — The uncertainty in redshift associated with the template, using the ±95% (full- 
drawn line) and ±90% (dotted line) confidence limits (the ±95% confidence limit is shown 
in Fig. 8). While Az/z is large for small redshifts, it drops to ~0.3 for z > 1.5 and even 
further at larger redshifts. Also shown is the Az/z coverage of 4 CO rotational lines with the 
planned ALMA receivers in the 3mm atmospheric windows. The redshift coverage depends 
on the observed frequency (the bandwidth being fixed at 16 GHz) and increases as a given 
transition is observed at higher redshift (lower observed frequency). Photometric redshifts 
are accurate enough to allow direct tuning for transitions at approximately z > 3. 



